Accepted to ApJ 

Preprint typeset using L^T^ style emulateapj v. 1 1/12/01 



COOL WHITE DWARFS REVISITED - NEW SPECTROSCOPY AND PHOTOMETRY 

Samir Salim'-^, R. Michael Rich'-^, Brad M. Hansen', L. V. E. Koopmans^, Ben R. 
Oppenheimer"* and Roger D. Blandford^ 

Accepted to ApJ 

ABSTRACT 

In this paper we present new and improved data on 38 cool white dwarfs identified by Oppenheimer et al. 
(2001) (OHDHS) as candidate dark halo objects. Using the high-resolution spectra obtained with LRIS on Keck 
I, we measure precise radial velocities for 13 white dwarfs that show an Ha absorption line. We show that the 
knowledge of radial velocities on average decreases the f/V-plane velocities by only 6%. In two cases the radial 
velocities put original halo candidates below the OHDHS velocity cut. The radial velocity sample has a velocity 
dispersion in the direction perpendicular to the Galactic plane of aw = 59kms"' - in between the values typically 
associated with the thick disk and the stellar halo populations. We also see indications for the presence of two 
populations by analyzing the velocities in the UV plane. In addition, we present CCD photometry for half of 
the sample, and with it recalibrate the photographic photometry of the remaining white dwarfs. Using the new 
photometry in standard bands, and by applying the appropriate color-magnitude relations for hydrogen and helium 
atmospheres, we obtain new distance estimates. By recalibrating the distances of the white dwarfs that were not 
originally selected as halo candidates, we obtain 13 new candidates (and lose 2 original ones). On average, new 
distances produce velocities in the UV plane that are larger by 10%, with already fast objects gaining more. Using 
the new data, while applying the same t/V-velocity cut (94kms"') and methods of analysis as in OHDHS, we 
find a density of cool white dwarfs of 1.7 x 10"^ pc"^, confirming the value of OHDHS. In addition, we derive 
the density as a function of the f/V-velocity cutoff. The density (corrected for losses due to higher f/V-velocity 
cuts) starts to flatten out at 150kms~' (0.4 x 10"^ pc""*), and is minimized (thus minimizing a possible non-halo 
contamination) at 190kms~' (0.3 x lO'^pc"-'). These densities are in a rough agreement with the estimates for the 
stellar halo white dwarfs, corresponding to a factor of 1.9 and 1.4 higher values. 

Subject headings: white dwarfs — stars: kinematics — Galaxy: halo — dark matter 



1. INTRODUCTION 

Dynamical studies indicate that the large quantities of gravi- 
tating matter exist in the haloes of galaxies, including our own. 
The amount of this matter far surpasses the matter with de- 
tectable electromagnetic radiation. Revealing the nature of this 
so called "dark matter" is one of the crucial goals in present-day 
astronomy and cosmology. 

One of the first proposed methods to indirectly detect dark 
matter was to look for lensing of background stars by the dark 
objects in the halo (Paczynski 1986). This technique, mi- 
crolensing, is sensitive to macroscopic objects, known as MA- 
CHOs, ranging in mass from planetary to stellar. At the time 
when the first microlensing experiments began, there was still 
no consensus as to whether the dark matter haloes should be 
baryonic or non-baryonic. With the exception of primordial 
black holes, microlensing would exclusively detect the bary- 
onic matter. Now, after years of observing, microlensing has 
ruled out a dark halo made entirely of MACHOs, but has 
nonetheless found more lensing events than expected from the 
known stellar populations, either in the Galactic disk and halo, 
or in the LMC (or SMC), where the lensed stars reside. The 
latest estimate from the MACHO group favors a halo in which 
dark objects comprise 20%, each with a typical mass of 0.6Mq 
(Alcocket al. 2000). Another microlensing experiment, EROS, 
puts upper limits for compact halo objects with stellar masses 



of 30%, from monitoring the LMC (Lasserre et al. 2000), and 
25% from SMC (Afonso et al. 2003). Some researchers (e.g., 
Sahu 1994) believe that lenses responsible for these events re- 
side in the LMC (or SMC) itself, however these scenarios have 
their own problems. 

Since the search for dark matter in the form of MACHOs be- 
gan, our paradigms about what the cosmological dark matter 
should be have changed. With the discovery of cosmological 
acceleration due to "dark energy" (Garnavich et al. 1998; Riess 
et al. 1998), the discrepancy between the critical mass density 
required by the inflationary model, and the low observed den- 
sities of gravitating matter was reconciled. Cosmological mod- 
els preferred the genuine dark matter to be in the form of non- 
baryonic cold dark matter, concentrated in galactic haloes. On 
the other hand, the large fraction, or all of the baryonic matter 
that was previously unaccounted for, is now believed to reside 
in the systems such as those responsible for the population of 
Ly a absorbers. 

This alone should be enough to shift the focus on the de- 
tection of non-baryonic halo dark matter (in the form of parti- 
cles), were it not for the MACHO results suggesting unknown, 
dark objects. For the problem to be more severe, there are a 
number of theoretical and observational arguments suggesting 
that the baryonic matter should at best constitute only a neg- 
ligible fraction of the dark halo, much lower than required to 



Department of Physics and Astronomy, University of California at Los Angeles, Los Angeles, CA 90095, samir@astro.ucla.edu, rmr@astro.ucla.edu, 
hansen@astro.ucla.edu 

^ Visiting Astronomer, Lick Observatory 

' Space Telescope Science Institute, Baltimore, MD 21218, koopmans@stsci.edu 

■* Astrophysics Department, American Museum of Natural History, New York, NY 10024, bro@amnh.org 
^ Theoretical Astrophysics, California Institute of Technology, Pasadena, CA 91125, rdb@tapk.caltech.edu 



2 



explain the fraction of it believed to be in MACHOs (Freese, 
Fields, & Graff 2000). Of the various candidate counterparts to 
MACHOs, ancient, cool white dwarfs are nevertheless the best 
ones, especially since their expected masses are comparable. 
The main question then, is, whether we can detect these halo 
white dwarfs, and if so, what fraction of the dark halo mass 
do they constitute. Ideally, these white dwarfs should account 
for all MACHOs; otherwise yet another form of baryonic dark 
matter is required. 

Surveying 10% of the sky, and using a proper- motion se- 
lected sample of objects with high reduced proper motions, Op- 
penheimer et al. (2001) (hereafter OHDHS) have identified and 
spectroscopically confirmed 98 white dwarfs, of which 38 have 
halo-like kinematics based on the derived velocity components 
in the plane of the Galaxy (UV plane). The derived mass den- 
sity of these cool white dwarfs indicated that they make up at 
least ~ 2% of the local dark-matter halo density, an order of 
magnitude higher than expected from the population of the stel- 
lar halo (as opposed to dark halo) white dwarfs (WDs). Stellar 
halo WDs should differ from the dark halo ones in their origin, 
since the dark halo white dwarfs were presumably produced 
in an early burst of star formation, and were cooling ever since. 
The OHDHS paper and its results stirred the astronomical com- 
munity, generating numerous objections - regarding the selec- 
tion of objects, kinematical cuts, contamination from non-halo 
populations, and so on. A thorough review on the subject of 
cool white dwarfs, and of various interpretations of OHDHS 
work is given by Hansen & Liebert (2003). 

In this paper, we re-analyze the OHDHS sample of cool 
white dwarfs using newly acquired high resolution spectra and 
CCD imaging. The new data allow measurement of radial ve- 
locities and determination of photometric distances in a more 
direct and precise manner, thus addressing some issues raised 
with respect to original OHDHS data. (The OHDHS used 
classification-grade, low resolution spectra, and photometry de- 
rived from photographic plates). 

In §2 we present spectroscopic observations, and the deriva- 
tion of radial velocities for white dwarfs exhibiting spectral fea- 
tures. Photometry from the CCD imaging is presented in §3, 
together with a calibration of photographic magnitudes used by 
OHDHS. In §4 we assemble all data to produce a new dataset 
required for kinematical analysis. Finally, in §5 we discuss the 
properties of the revised kinematical dataset, and use it to ad- 
dress some of the issues raised against the original OHDHS 
sample and its interpretation. In particular, we recalculate the 
densities of WDs using the OHDHS velocity cut, but also as a 
function of the cutoff velocity. In a forthcoming work the new 
dataset will be analyzed with additional techniques, including 
the kinematical modeling of stellar populations. 

2. SPECTROSCOPY 
2.1. Observations and data reduction 

The spectra of OHDHS cool white dwarfs were taken with 
the LRIS (Oke et al. 1995) spectrograph on Keck I, on three 
nights: 2002 Sept 13, 14 and 2002 Dec 4 (UT). The primary 
goal was to obtain high resolution spectra in the region around 
the Ha line (A6563). Thus, the spectroscopic setup in the red 
arm of LRIS consisted of 1200/7500 grating, giving disper- 
sion of 0.63 A pix~', and a spectral coverage of 5850-7140 
A. The only exceptions to this setup were in the cases of WDs 
* We use the same white dwarf designations as in OHDHS 



with peculiar spectra (that happen to show no hydrogen fines): 
WD2356-209^ and LHS 1402. The first was imaged in both 
the hi-res and the low-res modes, while the second only in the 
low-res: 300/5000 grating, 2.46 A pix~' dispersion, and 5010- 
10030 A range. The slit width of 1" gave the effective reso- 
lution in hi-res mode of 2.9 A. On 2002 Dec 4 one additional 
WD was observed (WD2346-478), this time with the 831/8200 
grating (0.92 A pix~' dispersion, 5630-7500 A range). 

The cumulative exposure times varied from 6 to 45 min- 
utes. Calibration lamp spectra were obtained at each pointing, 
and the internal-lamp flat-field images were taken once a night. 
Standard extraction and calibration IRAF tasks were employed 
to produce the final spectra. 

2.2. Measuring radial velocities 

Measuring the radial velocities was the main goal of the hi- 
res spectroscopy observing program. Targets were selected by 
inspecting the low-res spectra obtained by OHDHS. They found 
an Haline in 14 out of 38 cool WDs (denoted with an asterisk 
in OHDHS Table 1). Of these 14, all but three can be observed 
from Keek's latitude, and we have obtained spectra of all 11. 
In addition, we took spectra of another 12 WDs, thought to be 
featureless. Among these we find an additional two with an 
Ha line. Thus our radial velocity sample comprises of 13 cool 
WDs. 

The wavelength region covered by the hi-res spectra would 
allow detection of lines other than that of hydrogen (such as He 
and C). However, careful inspection of the spectra (with S/N 
ranging from ~ 30 to ~ 1 10) does not reveal any such lines. 

Each individual spectrum was wavelength calibrated against 
a lamp spectrum. The typical RMS of the calibration was 0.02 
A. The wavelengths of lines were first measured on individual 
spectra (of a same object) in order to evaluate the stability of the 
zero point of wavelength calibration, and to establish the accu- 
racy with which the central wavelength of Ha could be deter- 
mined. All measurements were done using the "splot" routine 
in IRAF. White dwarf spectra were normalized by the contin- 
uum, and a Lorentzian function was used to fit the lines. 

The stability of the zero point of the calibration was de- 
termined by measuring a bright sky emission line of [OI] at 
6300.304 A. The measured mean wavelength was 6300.30 ± 
0.01 A, indicating no systematic shifts in the calibration, while 
the scatter around the mean of 0.07 A (equivalent to 3kms"') 
gives the level of radial velocity error due to the wavelength 
calibration. 

The spectra of white dwarfs have different levels of signal to 
noise ratio, leading to variations in the quality of the Ha profile. 
In well exposed spectra, the non-LTE core of Ha line will be 
well defined, so measuring the central wavelength of the core 
would be superior to fitting a profile to the entire line. However, 
in lower S/N spectra the core may be degraded by noise. In or- 
der to test which method is more appropriate for our sample, we 
compared line fitting to the entire profile (6400-6800 A) (which 
was fitting the wing portions of the profile well, but often failed 
to fit the core), to the fitting of the non-LTE core alone (width 
~ 5 A). For each method, we first find the difference of the in- 
dividual measurements with respect to the mean value (for the 
given object), and then calculate the overall scatter. For wide- 
profile fitting the scatter is 0.45 A, while for core fitting it is 
0.32 A. In other words, core fitting seems to have better re- 
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peatability, and should thus be more precise. We should also 
note that we do not observe cases of split cores or of emission 
in the cores. 

Since the scatter, i.e., the error of the individual measure- 
ments, is much larger than the stability of the zero point of 
wavelength calibration, we conclude that it is safe to combine 
the spectra belonging to the same object (between three and 
six), thus obtaining a higher signal and eUminating deviant 
points by performing sigma clipping. Since the spectra of a 
given object were taken over a short period of time, the helio- 
centric velocity correction need not be applied at this stage. We 
then proceed by finding central wavelengths of Ha lines in the 
combined spectra, again by fitting a Lorentzian to the core. This 
gives our final measured wavelength. In order to evaluate the 
measurement error, we evaluate, as a function of a measured 
flux, the RMS scatter of the central wavelengths of individual 
spectra belonging to a given object. Not surprisingly, the scatter 
is larger for objects whose individual spectra had low intensi- 
ties. We find a linear relation between the logarithm of flux 
and the RMS scatter of individual wavelength measurements. 
However, at a certain flux level the RMS reaches a minimum 
value of 0.13 A, despite the increase of signal. In the fluxes of 
the combined spectra this plateau is actually reached for most 
spectra in our sample. Thus for most spectra, the total error 
from noise and wavelength calibration uncertainties is equiv- 
alent to w 7kms~\ well within the Umits acceptable for this 
study. Finally, to get the measured radial velocity, we apply 
the heliocentric correction. The observed radial velocities of 13 
WDs and their errors are listed in Table I. Also hsted are 10 
WDs observed with LRIS, but lacking spectral features. 

3. PHOTOMETRY 

3.1. Observations and data reduction 

Photometry was performed on CCD images taken at the 1 
m Nickel Telescope at the Lick Observatory, on 2002 Nov 27, 
Dec 3 and 4 (UT). "Dewar#2" CCD wifli a high quantum ef- 
ficiency extending to blue wavelengths was used. The first 
two nights were fully photometric and many standards over a 
large range of airmasses and colors were observed. This al- 
lowed construction of photometric transformations with linear 
and quadratic color-terms. The photometric accuracy from cal- 
ibration is 0.01-0.02 mag in aU bands. 

Of the 38 cool WDs, 19 reach an altitude high enough to be 
observed with this telescope. Of them, 18 were imaged in V 
and / (Cousins) bands. In addition, 9 of those were also ob- 
served in B band, and further 3 in ^ band. Of the 18 WDs with 
VI photometry, 9 also belong to the radial velocity sample. 

Object photometry was performed with an aperture equal 
to 1 FWHM of file PSF (typically I'.'A), and then aperture- 
corrected using a bright isolated star in the field. All /-band 
images were corrected for fringing. For each measurement a 
photon error from the object was combined with a photon er- 
ror of the aperture-correction star. The individual measure- 
ments were transformed into standard magnitudes and com- 
bined (weighted by photometric errors) into a single magnitude 
per star per band. Median photometry errors of the WD sample 
are (<7b) =0.053, (ctv) =0.035, (o-r) =0.030, (ct/) =0.035. The 
median error of V - / color, which we will use to deduce dis- 
tances, is 0.052 mag. The photometry is summarized in Table 
2. 

Since most of the cool WDs were not known before, it is not 
surprising that the literature search for photometric data pro- 



duced prior measurements for only two WDs: LHS 542 and 
LHS 147. The comparison of their photometry to ours is given 
in Table 3. It is in excellent agreement. Photometric magni- 
tudes in other bands do exist for several WDs in SDSS DRl, 
and for nine WDs in the 2MASS All-Sky Point Source Cata- 
log. The 2MASS measurements will be discussed in §4.3.1. 

3.2. Photometric calibration of OHDHS magnitudes 

Since we obtained CCD photometry for only one half of 
the OHDHS sample, it would be useful to derive photometry 
of other objects in standard bands. Thus, we would like to 
construct empirical transformations between the photographic 
plate magnitudes used by OHDHS: Bj, R^gv, and /n, and the 
standard photometric bands. Empirical transformations be- 
tween photographic and standard magnitudes do exist in the lit- 
erature, while the synthetic transformations can be constructed 
using the model spectra and the transmission curves, yet the 
first method has not been specifically appUed to stars such as 
cool WDs, while the second suffers from often ill-defined prop- 
erties of the actual response of a given plate/filter combination. 

Here we derive relations between photographic and standard 
magnitudes as measured by CCD photometry. 

ForSj: 

Bj=B-0.85(B-y) + 0.26. (1) 

This relation has a = 0.10. Since we know that as = 0.05, 
this indicates that gb, = 0.08, which is quite remarkable for the 
photographic photometry. Note that a high color-term indicates 
that (at least for WDs), Bj is actually closer to standard V than 
to standard B. 

For /?59f: 

/?59P = y-0.66(y-/) + 0.13. (2) 

Since we have our R for only 3 objects, we give this trans- 
formation relative to V. Excluded from the fit is LHS 1402 - 
a peculiar WD. Derived accuracy of the relation is 0.10 mag, 
while cTxjjp = 0.09. 

For/N: 

/N = /-0.09(y-/)-l-0.12. (3) 

This relation has a k, aj^ =0.16. Excluding objects for which 
OHDHS derive /n spectrophotometrically does not change the 
above relation. 

Another source of photographic magnitudes is the recently 
completed USNO-B catalog (a similar catalog, GSC-2.2, does 
not go deep enough in most cases). We have matched all the 
objects to counterparts in USNO-B catalog (Monet et al. 2003), 
and repeated the above analysis against B2, Ri and /serc - sec- 
ond generation sky survey magnitudes from USNO-B. How- 
ever, we find that USNO-B magnitudes are significantly inferior 
to those used by OHDHS, despite the fact that they come from 
similar or same plate material. Namely, we find (Tbj = 0.41, 
ctr^ = 0.58, and cr/sE^f, = 0.26. 

Overall, we conclude fliat flie OHDHS (that is, SuperCOS- 
MOS Sky Survey from which it is taken, Hambly, Irwin, & 
MacGillivray 2001; Hambly, Davenhall, Irwin, & MacGillivray 
2001) photographic plate photometry is of excellent quality, 
which lends credence to transforming them into standard mag- 
nitudes in order to derive photometric distances. 

Since we will be obtaining distances from V magnitude and 
V-I color, we want to directly transform OHDHS magnitudes 
and colors to these. We have seen that fij is quite close to V, so 
we use that magnitude to obtain the transformations 
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V=Bj-0.23(Bj-7n)-0.17, cr = 0.12, (4) 

y-/ = 0.62 (Bj-/n)- 0.04, £7 = 0.08. (5) 

These calibrations were derived omitting both peculiar- 
spectra WDs (LHS 1402 and WD2356-209). 

4. THE NEW DATASET 

4.1. Radial velocities 

4.1.1. Gravitational redshifts 

The observed radial velocities (Table 1) were extracted as 
explained in §2. However, they do not represent the true ra- 
dial velocities, since WDs exhibit substantial gravitational red- 
shift. The exact redshift depends on the mass and the radius of a 
WD, which we do not know for individual WDs in our sample. 
However, it is known that the range of these values is relatively 
small, so for our purposes it is sufficient to adopt a common 
value for the redshift. From Reid (1996) we find that the field 
WDs have an average redshift of 28.6kms~', with a spread of 
6.5kms~'. We subtract this value from the observed radial ve- 
locities, and add the scatter to the radial velocity measurement 
error. The final values are listed in Table 4. 

4. 1 .2. Common proper motion binaries 

In a case where a WD has a conmion proper motion compan- 
ion that is a main sequence star, one can obtain a measurement 
of a true radial velocity of a white dwarf (whether it contains 
spectral lines or not) by simply measuring the radial velocity of 
the main sequence component. Such a measurement circum- 
vents the gravitational redshift correction. To this end, we have 
carried out a search for companions in the USNO-B catalog, 
which lists proper motions based on multiple plates. Within the 
2' search radius we find no candidate companions with proper 
motions compatible to those of the white dwarfs. 

4.1.3. Selection effects 

The original OHDHS selection of cool WDs was based on 
the U and V components of the velocity. In the absence of ra- 
dial velocities, they were calculated by assuming W = 0. Since 
our goal is to characterize this population by obtaining the third 
component of the velocity from the radial velocities, we should 
try to evaluate whether the subsample for which radial veloci- 
ties are measured is representative of the population as a whole. 
Here we will restrict ourselves to a question of whether the sub- 
sample is representative kinematically - in terms of its U and V 
velocities, based on which the cool WD sample was selected in 
the first place. 

OHDHS selected their sample by requiring the WDs to have 
a velocity above some threshold in the t/V-plane. This thresh- 
old was chosen as 2cr UV velocity of the thick disk population 

UV = v/f/^ + (V+35kms-i)2 > 94kms~^ (6) 

One way of characterizing if the subsample of 13 WDs with 
radial velocities is representative, is to compare its average UV 
velocity with the typical UV velocities of randomly selected 
subsamples of 13 WDs out of the total 38. 

In order to obtain the distribution of UV velocities of ran- 
dom subsamples, we run a Monte Carlo simulation that draws 
13 out of 38 WDs numerous times, and for each drawing 

^ Available from http://research.amiih.org/users/bro 



calculates its average UV velocity. The average is taken in 
two ways: as a straight average, and as a weighted aver- 
age, where weights are the corresponding maximal volumes 
{{UV) = (E(f^^/Vmax))/(E(1/Vmax))) in which a WD could 
have been detected in the OHDHS survey. As explained in 
OHDHS, for each individual object, Vmax is set either by the 
magnitude limit of the survey, or by the lower proper motion 
cutoff, whichever is smaller. In Figure 1, the solid line repre- 
sents the distribution of unweighted averages. The unweighted 
average of the radial velocity sample is 195kms~', and is indi- 
cated by the arrow. We see that it is on the high side of the dis- 
tribution (thus somewhat favoring fast objects), but well within 
the spread of the distribution. The weighted UV velocity dis- 
tribution (dashed line) has two peaks, two lower being domi- 
nated by proper-motion limited objects, and the higher by the 
magnitude-limited ones. That we see two peaks might actually 
be indicative of the fact that the OHDHS sample is composed 
of more than one population. We see that the weighted average 
of our radial velocity sample lies right at the proper-motion lim- 
ited peak. If there really are two different population, this might 
mean that our radial velocity subsample is primarily represen- 
tative of one of these populations - the one with intrinsically 
lower velocities. This will be further discussed in §5.3. 

In any case, the radial velocity subsample does not seem to 
be extreme with respect to the whole sample in terms of its UV- 
plane kinematics. 

4.2. Proper motions 

Proper motions enter into the kinematical dataset since, to- 
gether with a distance, they determine two components of the 
physical velocity. OHDHS proper motions come from Super- 
COSMOS Sky Survey (Hambly et al. 2001). Since the OHDHS 
sample consists of relatively high proper motion stars, the aver- 
age fractional error is small (7% from listed values), and is thus 
not going to dominate in the velocity error, especially since the 
distances were originally derived from plate photometry. So, al- 
though not as important as other recalibrations, we nevertheless 
carry out a comparison of SuperCOSMOS proper motions of 
OHDHS WDs with those from the USNO-B catalog. USNO-B 
combines a large number of plates to arrive at a proper motion 
solution, the errors of which are found to be reliable (Gould 
2003). In this comparison we use data from B. Oppenheimer's 
online table^, since unlike the published version, it contains the 
individual components of the proper motion error, just like the 
USNO-B catalog. 

USNO-B contains proper motions for 36 OHDHS WDs. The 
median errors of SuperCOSMOS proper motions are 3.5 times 
larger than those of USNO-B. We find no systematic differences 
in two proper motion datasets. The reduced between the 
two datasets is 0.8, indicating good estimate of errors. (Super- 
COSMOS proper motions were recently also found by Digby 
et al. 2003 to agree with proper motions derived by combining 
SuperCOSMOS and SDSS positions). There are four cases in 
which either of the components is discrepant at > 2a level. In 
all of these cases the listed error of USNO-B proper motions is 
rather large, and also larger than the SuperCOSMOS listed er- 
ror. Visual inspection of DSSl and DSS2 images confirms the 
SuperCOSMOS value. Again, this is in Une with Gould (2003), 
who found that when USNO-B errors have large values, they 
are usually underestimated. Thus, except in these four cases, 
and one other in which USNO-B error is significantly larger 



than the SuperCOSMOS, in Table 4 we mostly list USNO-B 
values, with a flag indicating the source of proper motion. 

4.3. Distances 

4.3. 1 . Colors and atmospheric composition 

In principle, the multi-band photometry allows a determina- 
tion of the temperature of a white dwarf and its atmospheric 
composition. In practice, we are often limited by the range 
of photometric measurements and their precision. Neverthe- 
less, construction of color-color diagrams can be useful in some 
cases. 

From our CCD photometry we can place 9 OHDHS WDs 
onto a BVI diagram (Figure 2). The solid and the dashed tracks 
correspond to theoretical colors for ^ = 8 white dwarfs with 
pure hydrogen and pure helium atmospheres respectively, taken 
from Bergeron, Saumon, & Wesemael (1995) and Bergeron, 
Wesemael, & Beauchamp (1995). Tracks go from 12,000 K 
on the blue end, to 4000 K. Filled symbols correspond to WDs 
showing HaUne (therefore of the DA type). Judging from the 
position in the diagram, DA WDs are consistent with the hydro- 
gen track, as expected. Of the three non-DA WDs, LHS 542 is 
clearly not consistent with a hydrogen atmosphere. As shown 
by Bergeron, Leggett, & Ruiz (2001), a He atmosphere rep- 
resents a good fit to LHS 542 optical and infrared photometry. 
More interesting is another non-DA white dwarf, WD2356-209, 
an obvious outlier in the color-color diagram. OHDHS have al- 
ready shown its spectrum suggesting that it had "no analogs". 
Our LRIS spectra confirm this, and so does the photometry 
- we see excessively blue B-V color for an extremely red 
V -I = 1.97.** It has been suggested (I. N. Reid priv. comm.) 
that the heavy blanketing in the blue part of the spectrum is due 
to an extremely broad Na I doublet (which would make this 
WD a DZ type). Indeed, in our low-resolution spectra we see a 
well defined dip around 5893 A. Thus the subdued flux in the 
Na I region is consistent with V -I being boosted, and B-V 
becoming blue. Recently a similar WD was found with an ex- 
tremely wide Na I absorption line - SDSS J1330H-6435 (Harris 
et al. 2003). The S/N ratio of SDSS J1330H-6435 spectrum is 
too low to confirm the presence of other lines characteristic of 
DZ WDs. Even in our 45 min low-res exposure we cannot posi- 
tively identify the Ca II triplet. It is possibly absent because of a 
very low temperature. Finally, the two hottest WDs in our pho- 
tometry sample are also in this diagram, and judging from their 
position in it, they seem to have a temperature of ^ 11. 500 K. 
Their R magnitudes are also consistent with this temperature. 

We also look for flie OHDHS WDs in the 2MASS AU-Sky 
Point Source Catalog. Nine are catalogued. Since they are at 
the limits of 2MASS detection, the infrared photometry is rel- 
atively crude. Actually, one is not detected in H band, and ad- 
ditional five lack magnitudes. Of these, we have CCD pho- 
tometry for only 3 (one of which is LHS 542, discussed above). 
Therefore, for the remaining WDs we use a calibration between 
V and fij (Eqn. 4), whose accuracy is comparable to that of 
2MASS magnitudes. We then plot V -J against J-H va Fig- 
ure 3. Again, we can see that all WDs with Ha are compatible 
with model pure hydrogen atmospheres (Bergeron, Saumon, & 
Wesemael 1995). Three do not seem to be consistent with H 
colors, of which one is LHS 542. The oflier two (WD0205-053 
and J0014-3937) are even cooler, with temperature ~ 4200 K. 

Therefore, based on optical and infrared color information, 
we can conclude that 4 WDs in the sample probably have 

* Note that because Bj bandpass is actually close to V, WD2356-209 did not stand 
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He atmospheres: LHS 542, WD2326-272, J0014-3937 and 
WD0205-053. In addition, Bergeron (2003), based on OHDHS 
photometry alone, finds fliat WD0125-043 and LHS 1447 are 
better fit with a He model. 

Moreover, there are 3 additional WDs with V -I > 1.4, a 
red color that hydrogen WDs cannot attain. Finally, there 
are two WDs with peculiar spectra (the previously mentioned 
WD2356-209 and LHS 1402), which wiU be treated separately. 

This leaves 1 1 cool WDs which, based on the above argu- 
ments, could have either hydrogen or helium atmospheres. In 
order to use the appropriate color-magnitude relation, we would 
Uke to classify these remaining white dwarfs as well. Bergeron, 
Ruiz, & Leggett (1997) have shown that all hydrogen WDs with 
T > 5000 K (y-/ < 1 .0) should exhibit absorption Hues. There 
are four such WDs without H lines, to which we thus assign He 
atmospheres. We also see that of 6 WDs Ukely to have He at- 
mospheres based on color-color diagrams, all but LHS 1447 
(He WD candidate based on Bergeron, Leggett, & Ruiz 2001 
analysis of the original OHDHS photometry) have V -I >^ 0.8. 
Therefore, for the remaining 7 WDs (with 1.1 <V -I < 1.3), 
we also adopt He composition, but will allow for the possibil- 
ity of misclassification. Atmosphere assigrmients are given in 
Table 4. 

4.3.2. Color-magnitude relations 

Currently, only LHS 542 has a good trigonometric parallax 

measurement, with accuracy of 12% (but see 5.3). Therefore, 
we have to rely on photometric distances. It is to this end that 
we have acquired precise V -I CCD photometry. This was ac- 
complished for half of the sample. However, the CCD photom- 
etry also allowed the remaining WDs with only photographic 
photometry to be transformed into standard magnitudes (§3.2). 
Note that the calibrations between the plate and the standard 
magnitudes that appear in the literature were not measured (or 
modeled) for white dwarfs (e.g, Bessell 1986; Blair & Gilmore 
1982), so their use might introduce systematic offsets. On the 
other hand, our calibration is direct, so should be free of sys- 
tematic effects, and is accurate enough that it will not dominate 
in the final distance error. 

Next we need to decide what color-magnitude relation 
(CMR) to use. Here we have a choice between using an em- 
pirically measured relation (from WDs with trigonometric par- 
allaxes), or model relations. 

The most widely used empirical CMR is the one based on 
Bergeron, Leggett, & Ruiz (2001) multi-band photometry of 
152 WDs with measured parallaxes. The sample is assembled 
from heterogeneous sources and is dominated by disk WDs. A 
linear weighted fit to these WDs (of both DA and non-DA type), 
produces a relation mV—I 

Mv = (2.72 ± 0.09) (y-/) + (12.39 ± 0.07) (7) 

This relation has a reduced = 27 indicating that we are 
sampling a range of WD masses, or that some measurements 
are affected by multiplicity. 

OHDHS used Bergeron, Leggett, & Ruiz (2001) empirical 
data to construct their CMR, and then used model spectra with 
model bandpasses to convert standard magnitudes into plate 
magnitudes. 

In this paper, for deriving the distances, we will assume con- 
stant WD mass of 0.6Mq and thus use model cooling curves 

t in the original OHDHS color-color diagram (their Figure 4.) 



6 



for hydrogen and helium WDs. Surely, readers can use our 
photometry and a CMR of their choice to arrive at different dis- 
tance estimates. 

For H atmospheres we use Bergeron, Leggett, & Ruiz (2001) 
model CMR for WD mass of 0.6Mq. Cooling tracks for other 
masses are practically parallel, thus changing only the zero 
point of the relation. 

Mv = 3.42(y-/)+11.7, (y-/<1.3) (8) 

which is obviously steeper than the empirical relation given 
in Equation 7. 

For helium atmospheres, again using the Q.6Mq coohng 
curve of Bergeron, Leggett, & Ruiz (2001), the following linear 
fit is appropriate for the color range of interest 

Mv=2.38(y-7) + 12.7, (y-/>0.8). (9) 

The helium relation, on the other hand, is somewhat less 
steep than the empirical relation, at least for this red-color re- 
gion. In our cool WD sample there are two WDs with likely 
heUum atmospheres that are bluer than the above range. For 
them, we directly read values of My from the He cooling curve. 

As mentioned in the previous section, for 7 WDs we simply 
assume He composition based on color. If this assumption were 
not true, what would be the error due to using a wrong CMR? 
For V-I < 0.9, H and He cooling tracks almost coincide, so 
there would be almost no difference. For redder WDs, the error 
due to possible misclassification is 

aM^(class)= 1.10(y-/-l) + 0.13, (y-/>0.9). (10) 

The CMRs given above were constructed for fixed masses of 
0.6Mq. Various studies agree that this is a typical WD mass. 
However, the spread of masses seems to be less well known, 
and differs by as much as a factor of four (for a review see Sil- 
vestri et al. 2001). To be conservative, we will assume a value 
close to the higher estimates, ctm = 0.2M©. This mass range 
translates into absolute magnitude uncertainty of 

aMy (H) = 0.44; ctm, (He) = 0.50, (11) 
which we wiU use with the above color-magnitude relations. 

4.3.3. Peculiar white dwarfs 

As noted previously, OHDHS cool WD sample contains two 
objects with peculiar properties. It would therefore be inappro- 
priate to assign My to these objects based on y - / color alone. 

LHS 1402 has an extremely blue color ofV-I = -0.372, 
(CCD photometry). Taken at face value, this would indicate an 
extremely hot WD, with a temperature in excess of 100,000 K. 
Alternatively, this could be a very cool hydrogen WD where 
the blue color is the result of collision induced absorption by 
the H2 molecule (Saumon et al. 1994; Hansen 1998). In a 
pure hydrogen model this would indicate a temperature of only 
~ 2000 K and My of 18-19 (Saumon & Jacobson 1999). How- 
ever, Bergeron & Leggett (2002) have recently analyzed two 
somewhat redder WDs with similar spectra, and concluded that 
strong infrared suppression is better explained using a mixed 
H/He model, where He dominates. Thus, using a mixed model 
with N(H)/N(Ue) = 10-^ and g = 8, for LHS 1402 we obtain 
T = 3000 K and My = 16.4. Note that LHS 1402 is discussed 
in Bergeron (2003) as possibly having a pure hydrogen atmo- 
sphere, and thus being extremely faint and close. However, our 



LRIS spectra do not show a dip at 0.8/xm, suggesting that the 
mixed model is a better explanation. In either case, this could 
well be the coolest white dwarf known. 

Another WD with a peculiar spectral energy distribution, 
WD2356-209, has already been discussed in terms of its pho- 
tometry. One cannot use its very red y - / = 1 .972 to derive 
My from CMR. It seems likely that its temperature is in the 
3500^500 K range (Bergeron 2003), and thus we conserva- 
tively assign My = 16.5 ± 1 .0 to this object. 

4.3.4. Derived distances 

Finally, we use the appropriate CMRs to get absolute magni- 
tudes and thus the distance estimates from y- / color - both for 
the CCD photometry sample, and for the sample with photom- 
etry caUbrated from plates, using Equations 4 and 5. To get a 
total error in absolute magnitude, we add in quadrature the error 
in My due to the uncertainty in y - / color, the uncertainty due 
to a possible range of WD masses (Eqn. 11), and the uncertainty 
due to possible misclassification, where appropriate (Eqn. 10). 
The error in V magnitude is mostly negligible and is correlated 
with the V -I error, so we ignore it in calculating the distance 
error. The resulting distances and their errors are listed in Table 
4. 

For LHS 542 we thus obtain a distance of 3 1 ± 7 pc, in agree- 
ment with the trigonometric parallax distance of 3 1 . 1 ± 3 .6 pc. 
Mostly because of allowing for a large scatter in WD masses, 
our estimate of a typical distance error is 24%. If the spread 
in masses is actually two times smaller (O.IM©), the distance 
error would also be approximately cut in half. Our main goal, 
however, was to eliminate possible systematic errors that would 
affect the kinematics, and thus the interpretation of results. In 
the future, the parallaxes should provide a definitive check. 

Next, we compare the distances obtained here with those 
listed in OHDHS. Taken together (but omitting the two pe- 
culiar WDs), our new distances are 16% larger than OHDHS 
(13% if only WDs with CCD photometry are considered). This 
difference is not just an overall offset. Namely, for small dis- 
tances, the two distance estimates agree well, but the differ- 
ence increases farther out, reaching on average 0.55 mag in 
distance modulus for the farthest stars (i.e., new estimates are 
30% larger). Plotting the difference against y-/ shows a very 
similar trend - new distances are larger for the blue (almost 
exclusively hydrogen) white dwarfs. This is suggestive of the 
fact that the difference arises from our use of a model CMRs, 
which for hydrogen WDs (Eqn. 8) has a steeper slope than the 
empirical CMR (Eqn. 7), a type used by OHDHS. To test the 
significance of this effect, we recalculated all distances using 
the empirical CMR (Eqn. 7). The trend still exists, but is three 
times smaller. 

4.4. New cool white dwarf candidates 

The OHDHS cool white dwarf sample of 38 stars was se- 
lected based on their high velocities in the UV plane. Does 
our recalibration of distances make some white dwarfs exceed 
the velocity cutoff that they previously were not able to reach? 
To answer this, we look at 60 WDs that were identified by 
OHDHS, but had UV < 94kms"^ This list was not pubUshed 
in OHDHS paper but is available online, at the previously men- 
tioned URL. Of 60, for 39 WDs both Bj and /n magnitudes 
are available, so we use Equations 4 and 5 to get V and V -I. 
For the remaining 21 WDs we obtain V and y- / from Bj and 
/?59F, again calibrated using our CCD photometry. We then find 
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provisional distances according to hydrogen CMR (Eqn. 8) for 
all objects with V-I< 0.8, and the helium CMR (Eqn. 9) for 
red objects. Such assignment is conservative, so that even if an 
incorrect CMR is used, the WD will not be excluded. 

In this sample we notice two objects with the predicted dis- 
tances of only 7 pc. We identify one as LHS 69, a known nearby 
non-DA WD (Bergeron, Leggett, & Ruiz 2001). Its trigonomet- 
ric parallax gives a distance of 8.1 ± 0.3 pc. It has a standard 
V = 15.71 (we predict 15.69). The second is also a known WD 
with a mixed atmospheric composition (Bergeron et al. 1994), 
LHS 1 126, at a trigonometric distance of 9.9 ± 1 .0 pc, and hav- 
ing V = 14.50 (predicted 14.37). Both of these cases provide 
additional support that our calibration of plate photometry and 
the method for estimating distances are not biased. 

In this sample with new distances we find an additional 
13 WDs, previously falling below the velocity cutoff, to have 
UV > 94kms~'. Predicted UV velocities of these new cool 
WD candidates reach as high as 157kms"' in the case of LHS 
4041. Actually, this white dwarf, and one other of the thirteen 
(WD 0252-350) have actually already been proposed, based on 
all three components of the velocity, as candidate halo mem- 
bers (Pauli et al. 2003). Another six objects are also listed in 
the literature - of which five are spectroscopically confirmed 
white dwarfs. The full data for these new candidate halo WDs 
are given below the dividing line in Table 4.^ 

The effect goes in the opposite direction as well - some of 
the original OHDHS cool WDs no longer have estimated UV 
velocities in excess of 94kms~^ . These will be discussed in the 
§5.2. 

Table 4 provides all data or measurements (with errors when 
applicable) that are required for the kinematical analysis. Be- 
sides the original 38 OHDHS cool WDs, we list 13 WDs that 
qualify after the recalibration. For consistency with the original 
paper, the names of previously unnamed WDs are constructed 
from OHDHS J2000 coordinates. 

5. DISCUSSION 
5.1. Radial velocities and the UV -plane velocities 

OHDHS selected their sample based on the two components 
of velocity projected onto the sky. The third component, radial 
velocity, was not available. In order to obtain components of 
motion in the Galactic coordinate system {UVW), they assumed 
W = Q, which produces some arbitrary radial velocity that was 
then used to calculate U and V . Based on these U and V veloc- 
ities they selected their cool WD sample. Using the assumption 
oiW = Q, rather than Vrad = 0, was seen as a potential source of 
bias (Reid, Sahu, & Hawley 2001; Silvestri, Oswalt, & Hawley 
2002). The justification given by OHDHS is based on the fact 
that their sample is mostly in the direction of a Galactic pole, 
so U and V should not be much affected by the W component. 

For our radial velocity sample of OHDHS WDs we can cal- 
culate actual velocity components U, V and W, and thus di- 
rectly test and quantify the validity ofW = assumption. In Fig- 
ure 4 we plot as open squares the original"^ OHDHS positions 
of 13 WDs with radial velocities, and as filled squares their 
positions when the radial velocities are take into account. To 
help match the corresponding points we put sequential numbers 

' Here we note that H. Harris (2003, priv. comm.) pointed out that one of the low velocity objects, WDOl 17-044, is actually an almost equal WD binary, separated by 
3", and has measured for component A: V = 18.14±0.03, B-y = 0.90±0.05, y-/=1.00±0.03, and for B:V = I8.I7±0.03,B-y = 0.95±0.05,V-/ = 0.96±0.04. 
This would make the distance estimate, and thus the velocity, larger. However, using the USNO-B proper motion it still falls somewhat short of the 94kms~' cut. 

Note that here, and in the entire paper, for LP 651-74 (line 36 in the OHDHS Table 1) we use BJ-R59F = 0.72, and thus the OHDHS distance of 39 pc, in accordance 
with the online list. 



next to them. We see that the changes range from negligible to 
moderately high (point 12 = LHS 4033 moved by 81kms~^). 
On average, the points have moved by 29kms~^ What about 
the change in UV velocity? First we see that as a result of ra- 
dial velocity knowledge, two objects that were just outside of 
the 94kms"' cut have moved inward. These two cases actu- 
ally represent the largest changes (-48 and -58kms"^), while 
the average change is just -5kms~^ On average, each indi- 
vidual UV velocity is smaller by 6% (with a scatter of 21%). 
Thus it seeins that this is a modest effect, and that the use of 
W = Q assumption in the absence of radial velocities is appro- 
priate for this sample. Note that in this comparison we kept the 
sky-projected (tangential) velocities the same, i.e., we used the 
original OHDHS values. 

5.2. New distances and the UV -plane velocities 

Newly determined distances will, through modified sky- 
projected velocities, directly affect the derived values of U and 
V velocity components. We already saw the result of this in 
§4.4, where the new distances produced 13 new cool WD can- 
didates with potential halo kinematics. Also, in §4.3.4, we saw 
that there is a systematic trend affecting large distances more 
than the nearby ones. Since on average we expect high UV val- 
ues to belong to farther objects, we would expect that this trend 
will be reflected in the UV plane. The revised proper motions 
will also be responsible for some change, albeit a very sUght 
one. 

In Figure 5 we show the new (filled squares), and the orig- 
inal (open squares; equivalent to OHDHS Figure 3) f/V-plane 
positions of the 38 OHDHS WDs. In both cases the radial ve- 
locities are neglected, i.e., W = 0. To avoid clutter, individual 
points are not labeled, yet we notice that the new UV veloci- 
ties tend to be higher, especially for already high values. The 
analysis shows that on average points have moved by 46 kms~' , 
while the UV velocity on average increased by 23kms"' (with 
maximum change being -1-171 kms"'). Each individual UV ve- 
locity is on average larger by 10% (with a scatter of 30%). 

Thus, the net effect of radial velocities and new distance de- 
terminations is that the average UV velocities are somewhat 
higher than the original ones. 

5.3. Velocity component perpendicular to the Galactic plane 

For our radial velocity subsample of 13 WDs we can deter- 
mine true velocities in the direction perpendicular to the Galac- 
tic plane - the W component. We plot the W values in Figure 6 
as a function of the UV velocity. We have also added two WDs 
from the newly qualified cool WDs that have radial velocities 
measured by (Pauli et al. 2003) (open circles). Omitting the 
two WDs with UV ^ 94kms"', we find the W dispersion of the 
LRIS sample to be aw = 59kms-'. Only one WD (LHS 4033, 
V-I = 0.06) exceeds lOOkms"' - reaching = -153kms-i. 
Without it, we would have aw = 44kms"'. Actually, Dahn 
et al. (2004) are about to publish trigonometric parallax for 
this white dwarf. It turns out that LHS 4033 is very massive 
(M = 1 .25M0), and thus some three times closer than estimated 
based on O.6M0 CMR. This measurement also affects its sur- 
face gravity, producing much larger gravitational redshift than 
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we assume for our WDs. Altogether, if one was to use this in- 
formation a posteriori (which can be strongly argued against), 
we would get ayf = 45kms~^ Anyhow, these values are in be- 
tween the values usually derived for thick disk and spheroid 
(halo) populations, 35kms"' and 94kms"' , respectively (Chiba 
& Beers 2000). That the radial velocity sample probes mostly 
what appears to be a lower velocity (and younger?) population 
was already indicated in §4.1.3. Indeed, all but three of the ra- 
dial velocity WDs are proper-motion limited. While one can 
formally calculate au and <tv, they are meaningless because of 
the UV selection that was apphed to derive the sample in the 
first place. We defer a more thorough analysis of the kinemat- 
ics and population analysis to a forthcoming paper. 

5.4. Space density of white dwarfs 

The derived space density, and thus the mass density, of 
OHDHS cool white dwarfs was their key result, and the one 
that stirred most controversy, since it is considerably higher 
than the expected stellar (as opposed to dark matter) halo WD 
density. Using OHDHS original data for 38 cool WDs, and 
applying the same 1 /Vmax technique, with a limiting magni- 
tude of /?55f =19.8, one derives n = 1 .4 x 10"^ pc"^ (close to the 
value given by OHDHS of 1.8 x lO-^pc"^). Using a typical WD 
mass of O.6M0, this corresponds to p = 8.3 x lO~^M0pc"^, or 
some six times higher than the canonical value of stellar halo 
WD mass density of p,. = 1.3 x lO'^Mopc"^ (Gould, Flynn, & 
Bahcall 1998). (Note that this often quoted canonical value is 
somewhat of an educated guess, and not a real measurement (A. 
Gould, private communication 2003)). 

What estimate of density would we get with our updated 
data? In our analysis we include all cool WD candidates 
from Table 4, using the radial velocity data where available. 
For completeness, we append this Ust with the remaining low- 
velocity WDs identified by OHDHS (UV ^ 94kms"'). For 
them, we use listed SuperCOSMOS proper motions and do not 
check for the availability of radial velocities in the literature. 
Looking at the entire WD sample of 98 objects will allow us to 
characterize the density not just for a sample with a fixed UV 
cut, but rather as a function of the f/V -cutoff velocity. 

Of 47 WDs with UV > 94kms"' only 12 are magnitude- 
limited. All others, including all 51 WDs with UV < 94kms"^ 
are proper-motion hnoited, that is, the maximum distance at 
which they could be detected is determined by the proper- 
motion lower limit of the survey (330masyr"'). For the magni- 
tude limit we use ^59^; =19.8 transformed for each object into 
a corresponding V'™ using the Equation 2. 

If we restrict ourselves to the original UV-cut of 94kms~\ 
the 47 WDs that make this cut yield 

n94 = 1.72 x IQ-^pc"^ (12) 
which is similar or slightly higher than the original estimate 
based on 38 WDs and /^^^ = 19.8. For this sample we have 
(V /Vmax) =0.51, suggesting that R^^f^ has an appropriate value. 

How sensitive is this estimate on the choice of f/V-cutoff ve- 
locity? In Figure 7 we plot cumulative number density starting 
from the highest UV values (dotted line). However, since at 
higher UV values we Umit ourselves to yet smaller portion of 
the f/y-plane, we need to correct it. We do it by finding (at each 



UV) a fraction of objects with Chiba & Beers (2000) halo kine- 
matics that get excluded by different UV cuts. (The correction 
factor is normahzed to 1 for UV = 94kms~\ in order to make 
the results directly comparable to OHDHS who did not perform 
this correction. In any case, the actual correction at this velocity 
is very small). 

The corrected densities are shown with a solid line. The 
lower panel monitors (V/Vmax) at each point, and appears con- 
sistent with 0.5 for the entire range of interest. The vertical thin 
line denotes the 94kms"' limit. Inward of this limit we have a 
rise of density due to non-halo populations. Actually, we see 
that this rise begins inward of UV ~ 150kms~^ Note that in 
the corrected plot, the density of objects with halo kinematics 
should be independent of the UV cut. The corrected density at 
UV = 150kms-i is 

ni5o=0.42x 10-^pc-^ (13) 

corresponding to 1 .9 pc- The minimum value of the attained 
density is 

n„^ = 0.31xl0-4pc-^ (14) 

at UV ~ 190kms"' . This minimum value corresponds to 1 .4 pc 
- within the uncertainty of the stellar halo WD density. Be- 
yond UV ~ 250 km s"', the density estimate rises again. It ex- 
ceeds when UV > 400kms~^ However, in this range, the 
density is based on just 2-3 highest velocity objects with huge 
correction factors. 

6. CONCLUSIONS 

We obtain precise radial velocities for the majority of 
OHDHS WDs with HaUnes. This makes it possible to mea- 
sure more precisely the U and V components of the velocity, 
and also allows the W component to be derived. We show that 
the radial velocities do not affect significantly the way in which 
OHDHS selected their cool WD candidates. Our l^-velocity 
dispersion lies between the typical thick disk and halo values, 
an indication of a mixed sample. 

Also, our new CCD photometry, and the recalibration of 
OHDHS SuperCOSMOS plate photometry, allow for more ro- 
bust distance estimates. 

Finally, with the new dataset, and applying the same meth- 
ods of analysis as in OHDHS, we confirm the densities of cool 
white dwarfs that they derived. However, many times lower 
densities (consistent with the stellar halo) are found if one 
adopts higher J/V-cutoff velocities. This new set of data facili- 
tates a more sophisticated analysis, which we plan to present in 
a forthcoming paper. 

We thank Didier Saumon for valuable discussions, and the 
referee Hugh Harris for sharing some unpublished measure- 
ments. We are indebted to Nigel Hambly and Andrew Digby 
for their work in constructing the SuperCOSMOS sample. This 
publication makes use of VizieR and SIMBAD Catalogue Ser- 
vices of CDS in Strasbourg, France, and data products from the 
Two Micron All Sky Survey, which is a joint project of the Uni- 
versity of Massachusetts and the IPAC/Caltech, funded by the 
NASA and the NSF. 



Note that Dahn et al. (2004) result on LHS 4033 would only slightly affect n^^ (increasing it by 5%), and would leave ni5o and nmin unchanged. 
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Fig . 1 . — Kinematical selection effects of the radial velocity subsample. The two curves shows the straight (solid) and the weighted (dashed) average UV-velocity 
(UV = V J/^ + (V + 35 km s"' )^) of random subsamples of 13 white dwarfs (out of 38 OHDHS cool white dwarfs). Compared to these distributions are the actual 
(straight and weighted) C/V-velocity averages of our radial velocity subsample, shown with arrows. 
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Fig. 2. — B — V vs. V —I color-color diagram of OHDHS cool white dwarfc with CCD photometry. Filled circle.s are white dwarfs exhibiting Haline (DA type). 
The solid and the dashed tracks correspond to theoretical colors for g = S white dwarfs with pure hydrogen and pure helium atmospheres, respectively. The models 
cover 4000-12,000 K temperature range, with triangles marking every 1000 K. Labeled objects are discussed in the text. 
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Fig. 3. — V-J vs. J-H color-color diagram of OHDHS cool white dwarfs present in 2MASS AU-Sky Survey. See Figure 2 for legend. Tracks terminate at 
10,000 K. 
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Fig. 4. — The effect of radial velocities on the UV-plane kinematics. We plot the original OHDHS positions of our radial velocity sample (open squares), and 
their positions in the UV plane after the radial velocities have been take into account (filled squares). Numbers should help match the con'esponding points. 
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-400 -200 

V (km s-i) 

Fig. 5 . — The effect of recalibrated distances on the UV velocities. Original velocities of 38 OHDHS cool white dwarfs are shown as open squares (corresponding 
to their Figure 3), while values obtained with new distances (and proper motions) are shown as filled squares. 
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Fig. 6. — Component of motion perpendicular to the Galactic plane. W velocities of 13 OHDHS white dwarfs with radial velocities measured by us are shown 
as filled circles. Two open circles come from the "additional" cool white dwarf candidates, and were measured by Pauli et al. (2003). Vertical line represents the 
[/y = 94kms"' cut. 



16 



] 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 n — I 1 1 r 




500 400 300 200 100 

UV-velocity (km s~^) 

Fig. 7. — Number density of OHDHS white dwarfs as a function of UV-cut velocity. In the upper panel we show the corrected (solid line) and the uncorrected 
(dotted line) cumulative density of all 98 white dwarfs observed by OHDHS (summing from the high-velocity end). The corrected Une takes into account that at 
high UV values we sample a smaller part of the UV plane. Vertical line represents the UV = 94kms~' cut. In the lower panel we monitor average (V/Vmax) as we 
move towards the lower UV velocities. Note that beyond ~ 300kms~' we only have a couple of objects. 
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Table 1 

Observed WDs and their radial velocities." 



No.'' Name v*" (kms"') 



2 


WnOI 53-014 


-50 ±7 


5 


LHS 147 


—15 ±7 


7 




-28 ±7 


g 


T HS 4047 


—24 ±7 


9 






10 


TV UViiii / III 




12 


T HS 4033 


206 ±7 


13 


LP 586-51 


-22 ±7 


14 


WD2242-197 


-8±7 


15 


WD0205-053 




17 


WD0125-043 




18 


WD2346-478 


75 ±7 


20 


WD0300-044 


141 ±13 


21 


WD0123-278 




24 


LHS 1402 




25 


LHS 1274 


81±7 


27 


WD0044-284 




28 


WD2214-390 


50±7 


30 


LP 588-37 


-127 ±7 


32 


WD()045-061 




33 


WD0225-326 




35 


WDOl 17-268 




36 


LP 651-74 


74±7 



"In cases where no features were de- 
tected, radial velocity is omitted 

''Numbering follows Table 1 of 
OHDHS 



Table 2 



Johnson-Cousins CCD photometry 



No." 


Name 


B 


V 


R 


/ 


"obs'' 


2 


WD0153-014 




18.646 ±0.020 




18.415 ±0.034 


0202 


3 


LHS 542 


19.473 ±0.053 


18.221 ±0.023 


17.545 ±0.025 


16.926 ±0.028 


2321 


5 


LHS 147 


17.985 ±0.018 


17.615±0.011 




17.169±0.016 


3202 


6 


WD2326-272 


21. 027 ±0.165 


19.922 ±0.051 




19.095 ±0.051 


2202 


7 


WD0135-039 




19.644 ±0.048 




19.083 ±0.065 


0402 


9 


WD2356-209 


21.206±0.109 


20.850 ± 0.075 




18.878 ±0.033 


2404 


12 


LHS 4033 


17. 162 ±0.020 


16.992 ±0.017 


16.987 ±0.030 


16.936 ±0.036 


3222 


13 


LP 586-51 


18.318 ±0.039 


18. 185 ±0.022 


18.141 ±0.044 


18.096 ±0.039 


1212 


14 


WD2242-197 


20.504 ±0.098 


19.659 ±0.045 




18.861 ±0.074 


1201 


15 


WD0205-053 




18.898 ±0.145 




17.257 ±0.024 


0101 


17 


WD0125-043 




19.820 ±0.076 




18.911 ±0.054 


0202 


20 


WD0300-044 


20.782 ±0.266 


19.862 ±0.049 




19.301 ±0.051 


1404 


24 


LHS 1402 




18.050 ±0.014 




18.422 ±0.027 


0205 


27 


WD0044-284 




20.022 ±0.060 




18.713 ±0.039 


0402 


30 


LP 588-37 




18.496 ±0.024 




18.365 ±0.039 


0303 


32 


WD0045-061 




18.203 ±0.015 




17.219 ±0.019 


0302 


35 


WDOl 17-268 




19.057 ±0.049 




17.944 ±0.030 


0303 


36 


LP 651-74 


18.033 ±0.021 


17.342 ±0.011 




16.568 ±0.015 


2202 



"Numbering follows Table 1 of OHDHS 
''Number of observations for each band: BVRI 
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Table 3 

Comparison with the published photometry 



Name 


B 


V 


R 


/ 


Reference 


LHS 542 


19.23 


18.15 


17.53 


16.99 


Bergeron, Leggett, & Ruiz (2001) 




19.49 


18.25 






Liebert, Dahn, & Monet (1988) 




19.47 


18.22 


17.55 


16.93 


This work 


LHS 147 


17.97 


17.62 


17.38 


17.16 


Bergeron, Ruiz, & Leggett (1997) 




18.09 


17.66 






Liebert, Dahn, & Monet (1988) 




17.94 


17.57 






Eggen & Sandage (1967) 




17.99 


17.62 




17.17 


This work 
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Table 4 

The new dataset - relevant data on cool white dwarfs 



No. 




R.A. 


Dec. 




t^5 




^\t^6 ) 


Flap A^ 


y 


V — / 


FlatT p"^ 






Flap r*^ 


Vrad 








(deg) 


(deg) 


masyr 


masyr 


masyr"^ 


masyr ' 




mag 


mag 




pc 


pc 




kms ^ 


kms ^ 


1 


F351-50 


11.33178 


-33.49130 


1860 


-1486 


51 


10 


u 


19.37 


1.54 





40 


10 


He 






2 


WD0153-014 


28.46448 


-1.39468 


64 


-398 


4 


6 


u 


18.65 


0.23 


s 


170 


36 


H 


-79 


9 


3 


LHS 542 


349.78956 


-6.21383 


-618 


-1584 


1 


5 


u 


18.22 


1.29 


s 


31 


7 


He 






4 


WD0351-564 


57.78907 


-56.45198 


265 


-1052 


20 


19 


o 


20.96 


1.49 


o 


88 


22 


He 






5 


LHS 147 


27.03805 


-17.20401 


-120 


-1106 


7 


6 


u 


17.61 


0.45 


s 


75 


15 


H 


-44 


9 


6 


WD2326-272 


351.54458 


-27.24632 


576 


-104 


4 


9 


u 


19.92 


0.83 


s 


112 


27 


He 






7 


WD0135-039 


23.89029 


-3.95502 


454 


-186 


7 


6 


u 


19.64 


0.56 


s 


160 


38 


H 


-57 


9 


8 


LHS 4042 


358.64586 


-32.35540 


422 


-46 


1 


10 


u 


17.41 


0.18 





105 


25 


H 


-52 


9 


9 


WD2356-209 


359.18788 


-20.91370 


-329 


-211 


32 


20 


o 


20.85 


1.97 


s 


74 


34 


sp 






10 


WD0227-444 


36.87318 


-44.38573 


268 


-217 


12 


18 


o 


19.82 


1.06 





83 


22 


He? 






11 


JOO 14-3937 


3.44777 


-39.62331 


-226 


-714 


17 


2 


u 


18.70 


1.25 


o 


40 


10 


He 






12 


T HQ 




9 SQ|aA7 
— Z.O0v>+ / 


614 




1 n 


Q 

O 


TT 

u 


1 A QQ 


U.UD 


c 
o 




22 




1 7Q 
I/O 


Q 


13 


T P SSfi SI 
I-ir JO\J J I 


15 53001 


—0 S4QS^i 


342 


—122 




3 


u 


18 19 


09 




172 


37 




—51 


q 


14 


WD9947-1 97 


•^40 43478 


—19.67841 


346 


62 


4 


4 


u 


19.66 


0.80 




111 


27 


fl 


—36 




15 


W UUZU J - U J J 




— S 9QS1ft 


956 


400 






u 


18 90 


1 64 




29 


g 


He 






16 


wnni nn-fi4s 


1 S 9nQS7 


U'+ . '+Ovr+ ? 


516 


190 


5 





u 


17.58 


0.58 


o 


60 


14 








17 


WDOI 9S Od"^ 

W U\J i Z,J-\JH-J 


21 27431 




250 


—44 




4 




19 82 


91 


g 


98 


25 


He 






] 8 


VV UZ.J'-f-VJ '-V 1 o 


JlJU.J IZIJl 


—4.7 SSOf^O 


—270 


—454 




3 


u 


17 95 


83 


Q 


48 


1 1 


fl 


47 


10 


1 9 




49 nSzlQfS 




436 


322 


Q 




u 


18 43 


44 


Q 


106 


25 


He 






20 






—4 42355 


272 


—280 


17 


19 


Q 


19 86 


56 


g 


177 


41 




1 12 


15 


21 


VVUUiZO-Z /Q 


9n 7*^^74 






124 


J 


ID 


TT 

u 




1 90 


rv 

yj 


OU 


ZD 


tie; 






22 


WD2259-465 


3AA.11112 


-46.46632 


404 


-158 


4 


7 


u 


19.71 


1.26 





64 


20 


He? 






23 


WD0340-330 


55.03620 


-33.01671 


494 


-330 


3 


5 


u 


19.94 


1.17 





77 


22 


He? 






24 


T 140? 


35 13432 


—98 91 fS46 


492 


-30 


3 


2 


u 


18.05 


-0.37 




21 


5 


sp 






25 


LHS 1274 


24.80995 


-33.81756 


580 


-24 


3 


11 


u 


17.34 


0.49 





62 


15 




52 


9 


26 


WD0214-419 


33.56203 


-41.85251 


320 


-96 


18 


19 


o 


20.08 


0.98 





102 


25 


He 






27 


WD0044-284 


1 1 .00892 


-98 403 1 3 


-78 


-360 


13 


3 


u 


20.02 


1.31 


s 


69 


23 


He? 






28 


WD2214-390 


333.64480 


-38.98522 


1006 


-360 


2 


8 


u 


16.14 


0.67 


o 


27 


6 


H 


21 


9 


29 


WD2324-595 


351.04227 


-59.46895 


124 


-576 


5 


7 


u 


16.90 


0.14 


o 


88 


21 


H 






30 


LP 588-37 


25.58649 


-1.39757 


108 


-344 


1 


7 


u 


18.50 


0.13 


s 


186 


40 


H 


-155 


9 


31 


WD0345-362 


56.38631 


-36.18446 


142 


-588 


18 


67 


u 


20.40 


1.45 





71 


17 


He 






32 


WD0045-061 


11.27623 


-6.13876 


104 


-676 


3 


3 


u 


18.20 


0.98 


s 


43 


10 


He 






33 


WD0225-326 


36.36950 


-32.63163 


310 


160 


36 


4 


u 


18.61 


0.40 


o 


118 


28 


He 






34 


WD2348-548 


357.19527 


-54.76280 


364 


-96 


22 


32 


u 


19.21 


0.98 


o 


69 


17 


He 






35 


WDOI 17-268 


19.46521 


-26.81428 


476 


42 


3 


3 


u 


19.06 


1.11 


s 


55 


15 


He? 






36 


LP 651-74 


46.80880 


-7.24976 


-193 


-436 


1 1 


10 





17.34 


0.77 


s 


40 


8 


H 


45 


9 


37 


WD0135-546 


23.91108 


-54.59108 


660 


108 


17 


3 


u 


18.91 


1.13 





51 


14 


He? 






38 


WDOlOO-567 


15.17948 


-56.77684 


293 


293 


6 


8 


o 


17.44 


0.60 


o 


55 


13 


H 






Al 


WD2221-402 


535.46835 


-40.19267 


316 


-238 


17 





u 


19.81 


1.09 





80 


20 


He 






A2 


WD2342-225 


355.56890 


-22.45330 


312 


90 


6 


2 


u 


19.41 


0.82 





90 


22 


He 






A3 


WD0007-031 


1.77802 


-3.11857 


224 


-390 


3 


4 


u 


18.44 


0.79 





64 


15 


H? 






A4 


WD2236-168 


339.06532 


-16.79833 


318 


-60 


3 


6 


u 


18.48 


0.75 





69 


17 


He 






A5 


WD2234-408 


338.72467 


-40.75506 


287 


-249 


17 


14 


o 


17.72 


0.47 


o 


76 


18 


H 






A6 


LHS 1044" 


3.55330 


-13.18362 


-554 


-708 


2 


2 


u 


15.78 


0.68 


o 


22 


6 


H 






A7 


J0424-4551' 


65.99036 


-45.84513 


-100 


-532 


15 


1 


u 


16.85 


0.73 


o 


34 


8 


H 






A8 


LHS 3917^ 


348.82831 


-2.16120 


584 


192 


3 


1 


u 


16.48 


0.50 


o 


41 


10 


He 






A9 


LHS 4041 


358.57837 


-36.56524 


26 


-664 


1 


2 


u 


15.46 


-0.02 





59 


14 


H 


-27 


3 


AlO 


JL 193 


7.85905 


-44.63682 


342 


28 


13 


10 


u 


16.89 


0.1 1 


o 


91 


24 


He 






All 


LP 880-451" 


1.78131 


-31.22642 


336 


-124 


3 


6 


u 


16.47 


-0.12 





108 


26 


He 






A12 


LHS 1076' 


6.66973 


-55.41222 


-294 


-450 


5 


10 


u 


15.16 


0.23 





34 


9 


H 






A13 


WD0252-350' 


43.65459 


-34.83158 


44 


-328 


2 


3 


u 


15.79 


-0.05 





71 


17 


H 


86 


2 



Note. — Objects below the line are the new cool WD candidates, coming from OHDHS sample, but not Usted in OHDHS paper. Coordinates are given for epoch and equinox 2000. Radial 
velocities are corrected for gravitational redshift and come from this work, except for A9 and A13 (Pauli et al. 2003). 

"Astrometry source flag. U = USNO-Bl.O, O = OHDHS. 

"Photometry source flag. S = CCD photometry from this paper, O = Calibrated from OHDHS (SupeiCOSMOS) magnitudes. 

"Color-magnitude relation flag. H = Hydrogen CMR (DA WD), He = Helitmi CMR, H? = Hydrogen CMR used. Ha insecure. He? = Helium CMR used, but could be a non-DA hydrogen WD. sp 
= Special. 

''As discussed in 5.3, Dahn et al. (2004) are about to pubUsh trigonometric parallax of this WD, showing it to be 30 pc distant and bringing down the redshift-corrected radial velocity to 76kms"' 
(H. Harris, priv. comm.) 

'DA-type, V = 15.89, V-/ = 0.67, TTnig = 51.3 ± 3.8 mas (Bergeron, Leggett, & Ruiz 2001). 

'^DA9.5 (Scholz et al. 2000). 

«DZ7.5, Villanova WD Catalog (online); V = 16.31, V-I = 0.49, irug = 37.5 ± 5.9 mas (Bergeron, Leggett, & Ruiz 2001). 

'"DB3, Villanova WD Catalog (online). 

'DA5, V = 15.14, Villanova WD Catalog (online). 

^Name from Pauli et al. (2003). 



